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Abstract 

In this note we consider the two-dimensional risk model introduced in Avram et 
al. [2] with constant interest rate. We derive the integral-differential equations 
of the Laplace transforms, and asymptotic expressions for the finite time ruin 
probabilities with respect to the joint ruin times rinax('iti, 112) and Tmi 
respectively. 
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1. Introduction and Preliminaries 

Ruin theory for the univariate risk model has been studied extensively, see Asmussen 
[1], Rolski et al. [10] and many recent papers. Contrarily, there are only a few research 
on multivariate risk models. Chan et al. [5j studied the following two-dimensional risk 
model 
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where for fixed i = 1 or 2, {Xij,j — f,2,...} are i.i.d. claim size random variables, 
{Xij,j = f , 2, . . .} and {X2j,j — 1, 2, . . .} are independent, and both of them are also 
independent of the Poisson process N{t). 

Cai and Li i3i studied the multivariate risk model 



^ Uiit) \ 



V 



(1.1) 



where {{Xi 



,Xs_„),n > 1} is a sequence of i.i.d. non-negative random vectors, 
and independent of the Poisson process N(t). The model (jl.ip was further studied by 
Cai and Li in [4]. 

Yuen et al. [11] discussed the bivariate compound Poisson model 





/ Mi(i)+M(t) 



M2{t)+M{t) 
i=l 



I 



where Mi{t), A'hit) and M{t) are three independent Poisson processes, Xi{Yi) are i.i.d. 
claim size random variables, {Xi,i > 1} and {Yi,i > 1} are independent and they are 
independent of the three Poisson processes. 

Li et al. [7] discussed the bidimemsional perturbed risk model 



Uiit) 
U2{t) 




(T2B2(t) 



where N{t) is a Poisson process, {{Xij,X2j),j > 1} is a sequence of i.i.d. random 
vectors, {Bi{t), B2{t)) is a standard bidimensional Brownian motion, and the three 
processes are mutually independent. 

Avram et al. [2] studied the two-dimensional risk model below 



Ui{t) 
U2{t) 




Sit), 



(1.2) 



where S{t) is a Levy process with only upward jumps that represents the cumulative 
amount of claims up to time t, and the paper focuses on the classic Cramer-Lundberg 
model, i.e. S{t) is a compound Poisson process. 
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In this note, we discuss the above two-dimensional risk model (jl.2p with constant 
interest rate. About univariate ruin models with investment income, there have been 
a lot of research. Please refer to the recent survey paper Paulsen [8] and the references 
therein. 

Now we introduce our model. Let r be a nonnegative constant, which represents 
the interest rate. Then our model can be expressed as follows: 

U,{t) = e'-'u, + c, f e'^'-^'^dv - 5, f e''^'~'"US^,, i - 1, 2, (1.3) 

^0 "'0 

where Ui are the initial reserves, ci are the premium rates, and < (5i,52 < 1 with 

N(t) 

(5i + (52 = 1. 5t is taken to be a compound Poisson process, i.e. St = "'fcj ^ — 0, 

fc=i 

where N{t) is a Poisson process with intensity A > and {ak, fc > 1} is a sequence of 
i.i.d. random variables independent of N{t). Denote by F the distribution function, 
by / the probability density function of Ufc, respectively. Let Ok be the arrival time of 
the fc-th claim. Then we can rewrite (jl.3l) as 

N(t) 

U,{t) = e^'u, + ^(e'-* - 1) - 5, ^ e-^^'-^-Vfe, i = 1, 2. (1.4) 

fc=i 

For k = 1,2,..., denote by Tk the inter-time between the (fc — l)-th claim and the 
fc-th claim. Then {Tfc, fc > 1} is a sequence of i.i.d. random variables with exponential 
distribution with the parameter A, and 9k — X)i=i "^i- 
Define two joint ruin times by 

T„,i„(ui, U2) := inf{t > 0| min{(7i(t), (72(i)} < 0}, 
T,„ax(wi,M2) := inf{t > 0| max{{/i(t), {/2(t)} < 0}, 

and the corresponding ruin probabilities 

■i/'min(wi,M2) := P{T,„i„(ui , 1*2) < Oo}, 
■i/'max(wi,W2) := P {Tm2oS,U\ , U2) < Oo}. 

As in [2], we assume that ci/di > 02/62- Then if ui/Si > U2/S2, the above two joint 
ruin probabilities degenerate into one-dimensional ruin probabilities as follows: 

V'min(wi,M2) = 1p2{u2) P{'5t < OO S.t. U2{t) < 0}, 
V'max(wi,'"2) = -01 ("l) := P{^i < 0° S.t. Ui{t) < 0}. 
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Refer to [2] for the deduction. Throughout the rest of this note, we assume that 
ci/(5i > C1I82 and ui/(5i < U2l?>2- 

Remark 1.1. For each i, we know that 



Denote U{t) — {Ui{t),U2{t)), u = (wi,'«2) and = {civ — 6iSv, C2V — ^2*5*1)). Then we 
have 

U{t) = e^'u + ^ e'^(*-^)dZ„ = e"-' {^+ e-'^'dZ^ . (1.5) 
Differentiating both sides of (|1.5p relative to i, we obtain 

= rU{t)dt + dZt. (1.6) 
Integrating both sides of (|1.6p relative to t, we get 

U{t)^U{0) + r [ U{s)ds+ [ dZs. (1.7) 



By (|1.7p and the fact that {t,Zt) — {t,cit — SiS(t),C2t — 62S{t)) are a three-dimensional 
Levy process, following Protter 9, Theorem 32] , we know that U{t) is a two-dimensional 
homogeneous strong Markov process. 

The rest of this note is organized as follows. In Section 2, we show the integral- 
differential equations of the Laplace transforms of the joint ruin times T'inin(i*i, M2) 
and T'inax('"i, M2) respectively. In Section 3, we provide two asymptotic expressions for 
the finite time ruin probabilities with respect to the joint ruin time T,iiax(iti, M2) and 
7111^(^1,^2) respectively. 

2. Integral-differential equation 



In this section, we establish the integral-differential equations of the Laplace trans- 
forms of the joint ruin times T,„in(wi,W2) and rmax(ui, U2) respectively. 
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2.1. The result about Tjnin{ui, U2) 

In this subsection, we consider the joint ruin time Tniin(ui, U2). For convenience, we 
denote T'min(i*i, W2) by t(ui,M2). Its Laplace transform is defined by 

^minK,"2,s) :=^[e-^"("^'"=)] , fors>0. (2.1) 

Then 

< *i„i„(ui,M2,s) < 1. (2.2) 

Now we have the following result. 

Theorem 2.1. For ^ < ^ and s > 0, the function ^'min(-, s) satisfies the following 
integral- differential equation 

Cl \ 9*min , / C2 \ 9*min + « ,t, 

"2 H -7{ *min 



7' / du2 r 

H / *min(ui - <5iZ,U2 - 52Z,s)/(2:)dz = (2.3) 



with the boundary condition 



4',ni„(«i,^^ii,s) = £;[e-^^^(^"^)l , (2.4) 

where f{z) is the probability density function of Uk and T2 is the ruin time of risk 
process U2{t). Furthermore, ^'min *s the unique solution of i2.3\) - (K^ . 

Proof. Existence: For any > 0, by considering the occurrence time Ti of the first 
claim, we have 

_gjg-sr(«i,n2)] ^ i;[e-^^("i'"2),Ti > /i] +£;[e-^^("i'"^\ri < h]. (2.5) 

For any t > 0, denote by the information of the two-dimensional risk process 
{{Ui{s),U2{s)) : s > 0} up to time t, and by dt the shift operator of the sample 
path, i.e. (f?t(a;))s = cos+t for any sample path lo — {uJs,s > 0). By the properties of 
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conditional expectation and the strong Markov property, we have 
£;[g-sr(«i,«.)^ j,^ > /i] = £;[e-''^("i'"=)l{Ti>M] 

= / e-^''*,ni„ (e''''ui + -(e"-'' - 1), e'-'^ua + -(e"'* - 1), s ) Ae-^"du 
V r r / 

= e-(^+«)"vl,,„,„ (e'-'^wi + ^(e'-'' - 1), e^'^u^ + ^(e'^'^ - 1), ,s) . 

\ r r / 



(2.6) 



For the second item on the right side of (|2.5p . we have 

£;[e-«T(«i,«2)^y^ < h] 



-sr(ni,«2)^y^ < /i,cri < '- A LLjlI L 

e-'\Ti<h,ai> ^ . (2.7) 



Si 



By the strong Markov property, we have 
E 



, , e'-^iui + ^^(e'^^i - 1) e''^iw2 + ^(e''^i -1) 



Si 



<Ji < — F ^ — 

Si 02 



;*mi„ fe'^*ui + -(e''* - 1) - Siz, e'-'u2 + -(e''* - 1) - 52Z, s) f{z)dz. (2.8) 



On the other hand, if ci > 
and thus 



e'-^im + £L(e'-^i-l) e'-^iM2 + ^(e''^i-l) 



<5l 



52 



then t(wi, U2) = Ti, 



, , e'-'^iui + ^^(e'^'^i - 1) e'-^iw2 + ^(e'^'^i -1) 



"^\Ti<h,ai > 



Si 



\u, + ^(e'-'^i - 1) e^^^U2 + f{e^^' - 1) 



Si 



'■'ui + -£i(e'-*-l) i!'''u2 + — (<!'"'-l) 



7(z)d2 



(2.9) 
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By (1231)- (HH), we obtain 



*inin(wi,M2,s) 

\ r r J 



■!''*"l + ^(<!''*-l) «''*"2+^(<!'''-l) 

, /■ ^ 

-Atj+ / -St 



^0 



x*„.in (e'-'ui + ^(e''* - 1) - <5iz, e''*7/2 + -(e''* - 1) - '^2^, s) f{z)dz 
\ r r J 



h 



oo 



^^^''dt e-^*/(^)d^. (2.10) 

L A H 



By the definition of ^min(': 

know that if z > ^ X A 

then ^'mi„(e''*iti + f-(e'~* - 1) - 5iz, e/'ua + f (e'^* - 1) - (^iZ, s) = 1. By virtue of this 
fact and letting y := e^^ — 1, gi :— ui + ^ and q2 := M2 + we can rewrite p. 101) by 



*min(ui,W2,s) 

= e-(^+")''*„iin (ui + qiy, U2 + q2V, s) 
Jo Jo 

x*,„i„ fe'^*ui + ^(e''* - 1) - Siz, e'^'u2 + -(e'^* - 1) - <52Z, s] f{z)dz. (2.11) 
It's easy to check that y t if and only if /i J, 0. Hence by (|2.1ip . we have 



lini*min(Ml + qiy,U2 + q2y,s) = ^'min(ui,U2,s)- (2-12) 

By (j2.1ip . for any h > and y = e*"'' — 1, we have 

Q _ ^min(Ml + qiy, U2 + ^2^, s) - ^min(Ml, ^2, g) 

2/ 

^-{X+s)h _ 

H *min(til + giy, "2 + 92y, s) 




x^min f e''*Mi + -(e^* - 1) - ^iz, e^*U2 + -{e^' - 1) - (^2^, /(z)dz 



Z.-C. Hu and B. Jiang 



*min(Ml + qiy, U2 + q2y , s) - *mm(ui, "2, s) 



y 

g-(A+s)/i _ 2 



h 



'^mmiui + qiy, U2 + q2y, s) 



1 



x*min fe"*iii + -(e"* - 1) - S^z, e^*U2 + -(e''* - 1) - <52Z, s) (2.13) 
\ r r / 

By (|2.12p . letting y t 0, ft. 4- in the above formula and noticing that p.2p assures the 

interchange of limitation and integration, we obtain 

gi^:; K q2^ Wmin 



A 



OO 



*min(wi - 5iz, U2 - S2Z, s)f{z)dz = 0. (2-14) 

Replacing qi and q2 in (|2.14l) by wi + ^ and U2 + ^ respectively, we obtain the integral- 
differential equation. When ui/5i — U2/62, the joint ruin model degenerates into a 
univariate model, and then by the analysis in [2], we get the boundary condition. 

Uniqueness: By using similar arguments in Gerber j6j and noticing (|2.1ip . we 
define an operator T by 

Tg{ui,U2,s) = e"(^+'')''5 (ui + gi?/,U2 + 92?;, s) 

rh rOO 

Jo Jo 

xg fe''*ui + -(e"^* - 1) - Siz, e'^*W2 + -(e"* - 1) - S2Z, s) f{z)dz, 
\ r r / 

for any h > 0. It can be easily seen that ^'min is a fixed point of operator 7", as 

7~^min = ^min- Also, for two different functions gi and 52 we have for any h > and 

s > 0, 

|Tgi - T52I 

< e"(^+'')''|5i (ui + qiy,U2 + <722/, s) - .92 ("i + '7iy,'«2 + q2y,s) 
+ J Xe-^'dt J e""* 51 (e''*ui + -3-(e''* - 1) - ^iz, e''*U2 + ^(e'^* - 1) - <52Z, sj 

-.92 fe'-'ui + -(e--* - 1) - (5iz,e'-*U2 + -(e''* - 1) - <52Z,s) 
\ r r J 

<e-(^+^^'^\\g,-g2\\oo+( t \e-(^+^^'dt \ II51-52II00 



^^^^ II51-52II00, (2.15) 



Joint ruin probabilities with constant interest rate 



9 



where || • ||oo is the supremum norm over (mi,W2) G R^. Therefore, 7" is a contraction 
and by Banach's fixed point theorem and (|2.2I) . the solution of (|2.3p - (|2.4p is unique. 

Remark 2.1. One way to obtain the Laplace transform ^min of the joint ruin prob- 
ability TiTiin{ui,U2) is to solve the above integral-differential equation (|2.3p - (|2.4p nu- 
merically. A natural question arise: 

Can we give an analytical representation for the solution to the equation ll2.3\) - (K4\ l 
in some special cases such as exponential claim sizes ? 

Unfortunately, even in the case of exponential claim sizes, we have not found the 
way to solve the equation (|2.3p - (|2.4p . 

2.2. The result about TmaxCwi, U2) 

Define the Laplace transform of T'max(wi, U2) by 



where f{z) is the probability density function of Ufc and ti is the ruin time of risk 
process Ui{t). Furthermore, ^'max is the unique solution of i2.3\} - l[^.16]) . 

Proof. The proof is almost the same with that of Theorem l2.1[ and we need only to 
notice the following three things: 

(1) In this case, (|2.7|) becomes 



«'max(Mi,?/2,s) := E e"^^— t^^'"^) , for s > 0. 



Then we have the following result. 



Theorem 2.2. For < ^ and s > 0, the function ^max(',',s) satisfies the same 
integral-differential equation \2. S\] with the boundary condition 




E[e 



■^-Sr{m,U2)^rj.^ < h] 




1) e^T^U2 + f{e 



1 



+E e"\Ti <h,ai> 



V 




1)1 



where t{ui,U2) stands for rmax(ui,U2)- 



(2) r(wi,M2) = Ti, if (Ti > 



Si S2 



10 



Z.-C. Hu and B. Jiang 



(3) If . > + V then vl._(e-ui + f (e^ - 1) - 

Siz, e^'u2 + f (e'-* - 1) - 52Z, s) = 1. 

We omit the details. 

3. Asymptotics for finite time ruin probabilities 

In this section, we consider the finite time ruin probability associated with T'max(wi, U2) 
and T'min('Ui, W2). The original idea comes from [71 Section 4]. 

Define Xi{t) := e^^*Ui{t)/Si,i — 1,2. Then {Xi{t), X2{t)) has the same ruin times 
and probabilities with {Ui{t),U2{t)). Denote :=|^,Pi:=^,i = l,2. Then by (fTil) 
and our assumptions, we have 

N(t) 

X,{t) = +p,(l - e-''*) - « = 1,2, (3.1) 

k=l 

where pi > p2,xi < X2- 

For T > 0, define 4'niax(a;i, xs, T) P{7;„ax(i5ia;i, ,52X2) < T}. Then we have 

«'max(xi,a;2,r) = P{3t < T s.t. Xi{t) < and X2{t) < 0}. (3.2) 

Alternatively, we can also define '^jnin{xi,X2,T) := P{T^ij^{SiXi, 82X2) < T} and 
get 

'J'min(xi,a;2,r) = P{3t < T s.t. Xi{t) < or A2(<) < 0}. (3.3) 

In the following, we will provide asymptotic results on both ^E'max(2;i, a;2, 2^) and 
^I'min(a;i, a;2, T) under some condition. 

3.1. Asymptotic result about TmaxCwi? 1*2) 

Let T > 0, 77 G N, and {Vk,k = l,2,...,r7} be a sequence of i.i.d. random 
variables with the uniform distribution on (0, T]. Denote by {Vi , . . . , V*) the ordered 
statistic of (Vi, Vn). It's well known that conditioning on {N{t) = n}, the random 
vectors {9i, On) and {Vj', . . . ,V*) have the same distribution. Assume that {T4, k — 
l,2,...,ri} is independent of {(Jk,k > 1}. Define Ft{x) — P\e~'^^^ai < x]. Then we 
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have 



k=l 



= p|^e-'^^'=afe >x| =Ft™(x), (3.4) 
where Ft*"{x) stands for the n-muhiple convolution of Ft{x). 

Theorem 3.1. IJ cJk has a regularly varying tail with P{<7k > x} = L{x)/x'^, where 
L is continuous, slowly varying, hm L{x) — oo, and a > 0. Then for any T > 0, we 
have 

Before proving Theorem 13.11 we need one lemma. 

Lemma 3.1. Suppose that Uk satisfies the condition in Theorem \3.1\ Then Ft has a 
regularly varying tail. 

Proof. By the independence of Vi and cti , we have 
F^{x) = P{e~''^Vi > x} 



1 f^L{e'-yx), ^ S{x) 



-dy 



where S{x) = ^ ^j^rp^dy, which together with the assumption that L is continuous 
and lim L{x) = oo implies that 

x^oo 

lim S{x) = oo. (3.6) 
By the change of variable, we get that 

For any t > 0, by p.6p . (\3.7\i and the fact that L is a slowing varying function, we 
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obtain 



lim 3S = lim ^ = lim - ^ 



(e'-'J> ^l'^-^-' _ Lie-'-l'x) ~ ) L(e-'^'x) 



lim 



oo L{e'-'^x) 

= 1. 



lim 



1 - (e'^'^)" 



1 - (e"^^)" 
Hence Ft has a regularly varying tail. 

Proof of Theorem 13. li By Lemma [3.11 and [T, Proposition IX. 1.4], we know that 
Ft is a subexponential distribution. By (13. 2p and p.lSp . we have 

*max(a;i,a;2,T) 

(N(t) 

= P {Yl o-fe > Xi +p^{l - e-''*),i = 1,2, for some t < T 

I k=l 



( N(T) ^ 

oo 

= ^P{iV(r)=n} 

n=0 

I e-'-^Vfc > X, + p,(l - e---^), ^ - 1, 2 



k=l 



NiT) = 



(3.8) 



If Xi +Pi(l — e ''^) > X2 +^2(1 ^ e ''"^), then by p.lSp and the assumption that 
X2 ^2^1, we obtain 



fe=i 



N{T) = n 



-rT\ 



k=l 



NiT) = 



FT*"(xi+pi(l-e-'-^)) 



(3.9) 



and X2 +pi(l — e ''^) > xi +pi(l — e ''^) > X2 +P2(1 — e ''"^) > X2, which implies 
that 



Ft*"(x2 +Pi(l - e-"^^)) < Ft*"(xi + pi(l - e-"-^ )) < Ft*"(x2) 
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and thus 



(3.10) 



Ft™(x2) - Ft*"(x2) 

Sine 6 is a subexponential distribution, by T, Proposition IX. 1.5] and p.lOp . it holds 
that 



liminf^-*"^-^_±l^(^--^"»^l. 
By Fatou's Lemma, (13. lip and [1, Proposition IX. 1.7], we have 



(3.11) 



E P{NiT) = n}FT"\x, +pi{l - e-^')) 



lim inf 



ri=0 



\TFt{x2) 



liminfV P{iV(T) = n} 



Tl = 

oo 



fT*"(a;i+Pi(l-e-''^)) fT*"(a;2) 

^ nmmi ^ 

Ft*''{x2) ^T(a;2) 



-^^P{iV(T) = n}hminf:^I^ 
\1 ^ XI ^oo FT(a;2) 



Tl = 

OO 



\- J2 P{N{T) = n}r 



XT 



ri=0 



AT 



E[N(t)] ^ 1. (3.12) 

By (13^, ([53)1 and ((XT^ and under the condition that xi +Pi(l - 6"*^^) > a;2 +^2(1 - 
e"*"^), we have that 

hminf"^— ^'^^'^)>1. 
2:1^00 XTFt{x2) 

If XI +pi(l - e-'-^) < a;2 + P2(l - e"'''^), then 



P 



^ e-'^«'=afc >x,+p,{l- e-'-^), i - 1, 2 



fc=i 



iv(r) = n 



= p 



'^^'■cTfc > 2:2 +P2(1 - e 



lk=l 



NiT) = 



FT*'\x2+P2{l-e-^^)) 



(3.13) 



Since Ft is a subexponential distribution and X2 > xi, by [TJ Proposition IX. 1.5] we 
have 



j.^ PT*"(x2+P2(l-e-'-^)) 

a.™oo JV^(X2) 



(3.14) 
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Now By (I3.8p . (|3.13p and (|3.14l) . similar to the arguments in (|3.12p . we obtain that 
under the condition that xi + pi{l — e^^'^) < X2 +P2(1 — e^*"^) 

. r '^maxixi,X2,T) 

hmmi = > 1. 

Hence we always have 

liminf^Mfil^>l. (3.15) 

a;i->-oo \TFt{x2) 

On the other hand, by the assumption that X2 > xi, and p.lSp . we have 
*max(a;i,a;2,T) 

r N{t) 

= P < ^ e-'^'^^ak > X, +p,{l - e"''*),i = 1,2, for some t < T 

[ k=l 
( N(T) 

I k=l 



P{N{T) =n]plYl ^"'■'''^fc > ^2 

n=0 I fc=l 



7V(r) = n 



oo 



= ^F{iV(r)=n}FT™(a;2). 

ri=0 

By Fatou's Lemma, the above formula and 1, Proposition IX. 1.7], we have 

'^niax{xi,X2,T) 

hm sup = 

XTFt{x2) 



< lim sup ^ 

OO 



E P{N{T) = n}^^T™(x2) 



\TFt{x2) 



<^Y.Pi^^T) = n} lim sup 

-j^ OO 



AT 

= ^E[N{t)] = 1. (3.16) 
It follows from (|XT5l) and (|3J6l) that (jM]) holds. 

3.2. Asymptotic result about Tmin(Mi) 1*2) 

By Thcorcm l3.1l we can easily obtain the asymptotic result for '^minixi, X2, T), which 
is formulated as follow: 
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Theorem 3.2. // Ufe has a regularly varying tail with P{(7k > x} = L{x)/x°', where 
L is continuous, slowly varying, lim L{x) = oo, and a > 0. Then for any T > 0, we 

a;— >oo 

have 



^min{Xl,X2,T) 

lim = = 1. 

2;2>a:i-s-oo XTFt{xi) 

Proof. First, for i = 1,2, define 



(3.17) 



^P^{x^,T) = P{3t<T s.t. X,{t) < 0}, 

i.e. ipi{xi,T){i = 1,2) represents tlie ruin probability of Xi{t)(i = 1,2) within finite 
time T. 

Notice the fact that 

P{3t < T s.t. Xi{t) < and X2it) < 0} 

= P{3t < T s.t. Xi{t) < 0} + P{3t < T s.t. X2it) < 0} - P{3t < T s.t. Xi{t) < or X^it) < 0}. 
Then by and dSS]), we have 

*max(a;i,X2,r) = il;i{xi,T) +^'2(2:2, 7") - *min(a;i, a;2, 7")- (3.18) 
By Lemma lXn Ft is a subexponential distribution. Then by T, Proposition IX. 1.5], 



for i = 1,2, we have 



hm ifell = i. 



X2>a;i->cx) XPFxiXi) 

By p.lSp . (|3.19p . p.Sp . and the fact that X2 > xi, we obtain that 

^min(a^i,a^2,T) 
XTMxi) 

MxuT) - XTF^{xi)+Mx2,T) - *max(a:i,a:2,T) 



(3.19) 



1 



\TFt{xi) 



< 



< 



Mxi,T) - XTFt{xi) 



\TFt{xi) 



^2{x2,T) - *max(a;i,a;2,r) 



XTFrixi) 



XTFt{x2) 
ip2ix2,T) - *max(a;i,a;2,r) 



Ft(x2) 



XTFt{x2) 
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